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We present a renormalization procedure for Polyakov loops which explicitly implements the fact 
that the renormalization constant depends only on the ultraviolet cutoff. Using this we study 
the renormalized Polyakov loops in all representations upto the 27 of the gauge group SU(3). 
We find good evidence for Casimir scaling of the Polyakov loops and for approximate large-A'^ 
factorization. By studying many loops together, we are able to show that there is a matrix model 
with a single coupling which can describe the high temperature phase of QCD, although it is hard 
to construct explicitly. We present the first results for the non-vanishing renormalized octet loop 
in the thermodynamic limit below the SU(3) phase transition, and estimate the associated string 
breaking distance and the gluelump binding energy. By studying the connection of the direct 
renormalization procedure with a generalization of an earlier suggestion which goes by the name of 
the QQ renormalization procedure, we find that they are functionally equivalent. 
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I. INTRODUCTION 

The proof of confinement in QCD is now literally a 
million dollar question [1]. There are many ideas about 
the direction in which such a proof lies. No matter what 
these ideas are, once they are properly formulated, they 
are always open to test by lattice techniques. One of the 
long-lasting ideas has been to examine a toy model of 
QCD for large number of colors, 0. Because of the 
enhanced symmetry, many quantities become amenable 
to study in this limit. Computations of corrections upto 
sub-leading order, 1/N , have been made for quantities 
such as hadron masses and pion-nucleon scattering, with 
intriguing results. Recently, by adding supersymmetries 
to large- QCD, very simple toy models have been con- 
structed which are amenable to analytical treatment us- 
ing the AdS/CFT correspondence Q. Much excitement 
has been generated by the plethora of predictions of such 
toy models, and there have been exciting speculation 
about their applicability to QCD. 

At large A^ the dynamics of quarks is secondary to that 
of the gluons, being suppressed by power corrections in 
A^. Hence lattice tests of these ideas have been made in 
pure gauge, or quenched, QCD. There have been inves- 
tigations of the string tension and its scaling with N 4 1 , 
the nature of the phase transition with changing A^ 5|, 
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and tests of approximate scale-invariance of finite tem- 
perature QCD [6|. In this paper we investigate certain 
ideas about Polyakov loops and their behavior at finite 
temperature which have developed in recent years based 
on large A^ and the AdS/CFT correspondence. We inves- 
tigate matrix models which could be expected to describe 
the high temperature phase of pure gauge QCD. Our re- 
sults put very strong constraints on the kinds of matrix 
models which may provide a description of pure gauge 
QCD. 

The Polyakov loop is closely connected with confine- 
ment since it is the order parameter for the transition 
from a confined to a deconfined medium. Various mod- 
els based on Polyakov loops have been proposed to de- 
scribe the transition to a quark gluon plasma phase and 
its properties at zero as well as non-zero baryon de nsity 

0, i i [nl, [11, [11, [11, [11, IIE III III Hi US mf. 

Furthermore the connection of SU(3) theory to the large 
A^ limit, in the mean-field approximation, was discussed 
in [1, [^. For a test of the reliability and comparison 
of these Polyakov loop models to pure gauge theory and 
QCD with dynamical quarks, a detailed knowledge of the 
behavior of the Polyakov loop in the fundamental and 
higher representations in those theories is needed. 
The Polyakov loop needs to be renormalized, since it has 
divergent contributions from the ultraviolet. In Section 
mil we present a renormalization procedure which explic- 
itly incorporates properties expected of a good scheme. 
This direct renormalization technique is naturally appli- 
cable to Polyakov loop expectation values in all represen- 
tations of the color group. The multiplicative renormal- 
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izations in different representations are closely connected 
if the Polyakov loops satisfy a property called Casimir 
scaling. We present tests of Casimir scaling in Section 
IIVI This leads on, in Section |Vl to a test of large- fac- 
torization at = 3. We find good evidence for both in 
a temperature range not too close to Tc- 
Next, in Section IVll we examine whether the renormal- 
ized Polyakov loops are described in an effective matrix 
model. By examining the renormalization scheme depen- 
dence of these quantities, we find that a single parame- 
ter variation of matrix models describes the temperature 
dependence of the Polyakov loops in various representa- 
tions. We also show that the matrix model is unlikely 
to consist of a small set of terms, and therefore hard to 
construct explicitly from the phenomenology of Polyakov 
loops. 

In Section lVIII we consider the adjoint Polyakov loop cor- 
relations below Tc- We report the first measurement of 
the renormalized adjoint Polyakov loops in the thermody- 
namic limit at finite temperature in the confined phase 
of QCD. We find that aspects of the adjoint Polyakov 
loop correlations can be summarized in the physics of 
gluelumps, i.e., colorless states made of static adjoint 
sources and glue. 

In the appendices we consider a renormalizati on p roce- 
dure, the QQ procedure, earlier suggested in ^22]. We 
extend it to the renormalization of Polyakov loops in ar- 
bitrary representations, consider the relation between the 
direct and QQ renormalization procedures, and examine 
"color averaged" Polyakov loop correlators in various rep- 
resentations. 



II. DETAILS OF THE CALCULATIONS 

We have performed simulations for temperatures up 
to 24 Tc on Nl x Nr lattices with Nr = 4, 6, 8 and 
up to 32 in SU(3) pure gauge theory with the tree level 
Symanzik- improved gauge action [2^, |23|. We used a 
pseudo heatbath algorithm with FHKP updating in the 
SU(2) subgroups. Each heatbath update was followed 
by four overrelaxation steps. The statistics varies from 
1000 to 10000 of such sweeps after suitable thermaliza- 
tion. The physical scale has been set using the zero tem- 
perature string tension, a, (25| and a determination of the 
critical coupling for the deconfinement transition from 
[2^ . We have calculated the Polyakov loops in all repre- 
sentations up to £) = 27 using the operators defined in 
(jA12[) - (|A18p . The errors on the observables were deter- 
mined with the Jackknife method. 

Furthermore, we have reanalyzed configurations gener- 
ated using two flavors of staggered quarks with mass 
m/T = 0.4 on a 16'^ x 4 lattice at several temperatures 
above and below the transition temperature [23, . At 
each temperature we have used statistics of several thou- 
sands to calculate Polyakov loops up to I? = 27. 



III. RENORMALIZATION OF POLYAKOV 
LOOPS 

We define the thermal Wilson line, P{x) at spatial po- 
sition X as 

P{x)= n %.^),o, (1) 

1=0 

where U^g^i)^^ is the gauge link matrix in the time di- 
rection at the point x and Euclidean time i. U is a 
3x3 matrix which belongs to SU{2>). We define the 
local Polyakov loop as the trace of P{x), 

L(x) = TrP(a-), (2) 

where the trace is normalized to one. We denote the ex- 
pectation value of the Polyakov loop in the fundamental 
representation of SU(3) by 

L^^{^Y.L{S))- (3) 

X 

Polyakov loops in different representations, L^, are de- 
fined in Appendix |^ The subscript D indicates the di- 
mension of the color (irreducible) representation of the 
Polyakov loop, e. g., L3 for fundamental or Ls for ad- 
joint. Expectation values of Polyakov loops are ultravio- 
let divergent. We will use the superscripts 6 or r for bare 
and renormalized Polyakov loop respectively. 

A. Basic properties of renormalization 

It was pointed out by Polyakov [2^ that for smooth 
loops, ultraviolet divergences can be absorbed in the 
charge renormalization of gauge fields: 

U^{T)^{ZD{g^)Y^^HU9% (4) 

where 1{C) is the length of the contour and the coupling, 
= 6//3, on the right is the bare coupling. The quan- 
tity on the left is properly renormalized and depends on 
the renormalized coupling or, through this running cou- 
pling, on the temperature. Cusps and self intersections 
of loops give rise to logarithmic divergences which de- 
pend, e. g., on the angle of the cusps ^2^,3(^,^31]. Spatial 
averages of operators such as Tr[L(a;)]", which wind n 
times around the lattice, also need separate renormaliza- 
tion. Similarly, composite operators such as powers of 
Polyakov loops, including Polyakov loop susceptibilities, 
also require independent renormalization. 
We have written eq. ^ for an arbitrary representation 
D. The renormalization constants in different represen- 
tations, Zo, can be related to each other if both the bare 
and the renormalized loops satisfy the relation 
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called Casimir scaling. Here €2(0) = Tr^^A^A" is 
the quadratic Casimir operator in the representation D. 
When Casimir scaling holds, the quantities 



l/do 



(6) 



where dj) — C2(I?)/C2(3), are all equal. 

In keeping with the general form in eq. ^ we can write 



(7) 



where the renormalization constants Zu(g'^) should only 
depend on the bare coupling. Such a multiplicative renor- 
malization is expected to compensate entirely for the de- 
pendence of the bare loop on the cutoff, i.e., the bare cou- 
pling, so that the renormalized loop on the left is a func- 
tion only of the temperature. Furthermore, if Casimir 
scaling is found to hold, then all the Zd collapse to a sin- 
gle function Z^(g^), i.e. loops in all representations can 
be simultaneously renormalized. Note that there is one 
remaining ambiguity: the function ^3(3^) can be multi- 
plied by a single coupling independent constant without 
affecting the renormalization. Thus, a one parameter 
family of renormalization schemes for Polyakov loops is 
defined by eq. ([7]). 



B. Direct renormalization of the Polyakov loop 
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FIG. 1: Comparison of the renormalized Polyakov loop in 
the fundamental representation, Uj,{T), obtained with the two 
different renormalization procedures. 

In this subsection we present a complete renormal- 
ization procedure which implements eq. We call 
this the direct renormalization prescription. A simi- 
lar procedure was discussed earlier in [34|- Denote by 
L£)((7^, Nr) the Polyakov loop expectation value obtained 
after taking the thermodynamic limit at a temperature 
T = l/a{g'^)Nr, where a is the lattice spacing at a bare 
coupling g^ and N^- is the temporal extent of the lattice. 




FIG. 2; Comparison of the renormalization constants for fun- 
damental loops, Z^{g^), for the two different renormalization 
procedures, denotes the bare coupling and the solid line is 
the same as in fig. |4] 



We will describe the procedure for a fundamental loop 
first. 

First choose the value of L'^{Tref), at a reference tem- 
perature Tref- It is clear from eq. ([7]) that this choice is 
exactly equivalent to fixing the renormalization scheme. 
It is convenient, but not necessary, to take Tref to be 
the maximum temperature in the study: in our case 
Tref — ^'2Tc. We discuss our choice of scheme in Ap- 
pendix [C] This sets the first step of the iterative proce- 
dure starting at the initial temperature Ti = Tref- 
Next we need measurements at (at least) two different 
temporal extents, A^" and A^^, say, with > , both 
at the temperature Ti (we begin with i — I and set up an 
iteration). Therefore these measurements correspond to 
two different lattice cutoffs a{gf^)N^ = a{gfij)N^ with 
ai,a < where the subscripts are self-explanatory. We 
obtain two different renormalization constants, 



{Z.igljf^ Lligl^,N:^) = Ll{T,), 



Z3{glp)f' Ll{gl,,N,^) = Ll{T,). 



(8) 
(9) 



The third step is to advance the iteration. We do this 
by making a measurement of L'^{gfp,N^) on the lat- 
tice with temporal extent A^" at a temperature Tj+i = 
l/ai^/jNr ~ (N^ /N")Ti. Since the renormalization con- 
stant is already known from eq. ([9]), one obtains the value 
of L5(T,+i) 



l) = (Zsiglp)) 



(10) 



Since we have the value of the renormalized loop at a 
new temperature, we can now iterate the procedure from 
the second step on. The iteration gives the renormalized 
loops and the renormalization constants at a decreasing 
series of temperatures. 

Four points about the prescription are worth noting ex- 
plicitly. First, the procedure extends without change to 
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FIG. 3: The renormalized fundamental Polyakov loop in 
SU(3) pure gauge theory for two values of the temporal lattice 
extend A''t. The lines show the perturbative result (|lip and 
are explained in the text. The arrow represents the asymp- 
totic high temperature limit, = 1. 



any representation D. The test of Casimir scaling would 
be to assume that the bare loops in different represen- 
tations at Tref are related by eq. (O, and then check 
whether the renormalized loops at all T are related in 
the same way. We discuss this further in Section IIVI 
Second, in the confined phase of the pure gauge theory 
the bare Polyakov loop, in any representation with non- 
vanishing triality, vanishes in the thermodynamic limit; 
as a result the direct renormalization procedure can only 
be used above Tc for such representations. Third, a 
reverse iteration can always be performed by choosing 
T^+^ = l/ai^aN^ = {N^/N^)T, > T,. Finally, although 
we discussed the procedure for two values of N-r, it can 
be easily extended to a larger number of values for the 
temporal extent. 

The renormalized Polyakov loop in the fundamental rep- 
resentation obtained by the direct procedure described 
here is shown in Figure [TJ Also shown, for compari- 
son, are the results obtained from a completely differ- 
ent renormalization procedure [2^ based on a matching 
of the short distance behavior of heavy quark-antiquark 
free energies to the zero temperature potential (labeled as 
QQ-renormalization). Both these procedures allow a one 
parameter family of renormalization schemes, and the 
schemes have been chosen so that the value of L^\Tref) 
match. FigureOshows the results for the renormalization 
constant. These figures indicate the functional equiva- 
lence of the two renormalization procedures. 



C. Fundamental Polyakov loops 

We have extended previous measurements of the fun- 
damental Polyakov loop [12] to temperatures as high as 



FIG. 4: The renormalization constants Z3{g^) as a function 
of the bare coupling = &/ f3 calculated on lattices of size 
32=* X Nr with iV^ = 4 and 8. The line comes from a fit to 
(fT3)) as explained in the text. 
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TABLE I: The renormalized fundamental Polyakov loop 
Ll{T) obtained on lattices of size 32^ x N-r with A^^ = 4 
and 8. T /Tc denotes the temperature in units of the critical 
temperature. 



24Tc. The results for L'^{T) are shown in Figure [3] and 
listed in Table[Tl The corresponding renormalization con- 
stants are plotted in Figure [3] and listed in Table [TTl 
The direct renormalization procedure for the fundamen- 
tal Polyakov loop stops at corresponding to Tc on the 
lattice with the smallest Nt- Since the QQ procedure 
gives identical results upto this point, and can be con- 
tinued to larger , the tables contain results obtained 
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TABLE II: The renormalization constants for the fundamen- 
tal Polyakov loop, Z^lg^) obtained on lattices of size 32^ X Nr 
with AV = 4 and 8. = 6//3 denotes the bare coupling. 



using this procedure. 

The Polyakov loop for SU(N) pure gauge theory in HTL 
perturbation theory [sl] is 



Ld = 1 



3^, 

2 X 2 
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87r2 
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1 /27r27V^ 3 



(11) 



With an appropriate running coupling, g{T), this defines 
the renormalized Polyakov loop up to 0{g'^). We make 
the specific choice of the two-loop formula. 



g-2(r) = 2/3oln 



jjT_ 

2 



Po \ U 



MS 



, (12) 



with Po = ll/(167if) and 5i = 102/(167r2)2. and 
Tc/Ajjg = 1.14. [H, m, ill. These predictions are 
shown in Figure [3] for the choices /i = tt/2, tt and 2tt. 



Due to the phase transition, the Polyakov loop expecta- 
tion value vanishes below Tc and rises beyond ST^. It 
starts to decrease from about lOTc and approaches the 
asymptotic high temperature limit, L"!^ — 1 (indicated by 
the arrow in fig. [3]), from above, in qualitative agreement 
with weak coupling theory. The lattice measurements 
seem to fall a little slower than the HTL prediction, upto 
the highest temperature examined. Approximate quali- 
tative agreement with HTL perturbation theory, without 
exact quantitative agreement upto very high temperature 
has been seen in many other contexts in high ternpera- 
ture QCD, most notably for screening masses [stI. \3§^ at 
high temperatures. 

The comparison of the renormalization constants Z^{g^) 
for different Nr, shown in Figured demonstrates that 
depends only on the bare coupling, and not on tempera- 
ture. The solid line in Fig. [4] shows the result of a fit with 
a (two-loop) perturbation theory inspired Ansatz, 



^3(5^) = exp g 



1 



N 



(13) 



where Q^^^ and Q^^^ are expected to be independent 
of Nr if eq. (O is to be satisfied. Although we are in 
coupling range which is not small enough for the weak 
coupling expansion to be numerically accurate, the fit 
works surprisingly well. In fact the bare coupling be- 
comes significantly larger than unity before this ansatz 
begin to overestimate the actual values of Z3. From the 
best fit analysis with a fit range g^<1.2 we obtained the 
values Q(2) ^ 0.0591(21) and ^ 0.0605(54). In- 

terestingly, although our computations are done using a 
Symanzik improved action, the value of Q*-^^ agrees rea- 
sonably well with result from lattice perturbation the- 
ory, Q(2) ^ 0.057(2) [11], using the Wilson action. The 
renormalization scheme dependence of these results will 
be discussed in Appendix [Cl 



D. Polyakov loops in other representations 

The results for the measured bare Polyakov loops is 
shown in fig. [5] for the representations D — 3 to D — 15. 
The renormalization of Polyakov loops in each represen- 
tation, D, was obtained using the direct renormalization 
procedure. For each D, the starting point was taken at 
Tref = 12Tc. We fixed the scheme through the choice 



(14) 



Except for this assumption at a single temperature, the 
renormalization was performed independently at each D. 
Since the loops were measured at arbitrary temperature 
values, spline interpolations (solid lines) for the data sets 
were used in the renormalization iteration. The errors 
on renormalized Polyakov loops and renormalization con- 
stants were obtained through a jackknife analysis. The 
accumulation of errors during iteration, the exponential 
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FIG. 5: The (bare) Polyakov loops for different represen- 
tations D measured on 32'^ x 4 lattices. The solid lines are 
the splines used in our analysis to extract the renormalization 
constants in the direct renormalization procedure. 



FIG. 7: The Casimir-scaled bare Polyakov loops for different 
representations D measured on 32^ x 4 lattices. 
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FIG. 6: The renormalized Polyakov loops for different rep- 
resentations D obtained with the direct renormalization pro- 
cedure. Also shown are the results obtained from the QQ- 
method for fundamental and adjoint loops, labeled Dqq. 



in ([7]), and the larger statistical errors for higher rep- 
resentations lead to large errors in the renormalization 
procedure for representations higher than Z? = 8 as one 
approaches Tc. 

The results for L^jy{T) are shown in Figure [S) For com- 
parison we have also included in the figure the results 
for the fundamental and adjoint representation obtained 
with the QQ procedure (see Appendix IB|) . These agree 
within errors, demonstrating again that the two renor- 
malization procedures give equivalent results. 
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FIG. 8: The renormalization constants for different represen- 
tations D obtained with the direct renormalization procedure. 
Also shown are the results obtained from the QQ-method for 
fundamental and adjoint loops. 



IV. CASIMIR SCALING 

In [iOl it was shown that Casimir scaling is realized 
in perturbation theory (at least) up to two-loop order, 
0{g'^). This statement even holds for QCD with (mass- 
less) dynamica l q uarks as shown in lattice perturba- 
tion theory in [4l|. Moreover, lattice calculations at fi- 
nite temperature employing an effective action for the 
Polyakov loop in SU(3) have found Casimir scaling to be 
realized for the Polyakov loop as well [i^. Numerical 
calculations on the lattice at T = in SU(3) pure gauge 
theory show that Casimir scaling is realized also in the 
non-perturbative regime for distances smaller than the 
string breaking distance (43j . The very good agreement 
of the lattice data with the Casimir scaling hypothesis at 
non-perturbative distances in the vacuum has consider- 
able ramifications on models for non-perturbative QCD, 
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especially for the confinement mechanism [44|| . 
We have noted before that if Casimir scaling holds, then 
it holds for bare as well as renormalized loops. Since bare 
loops have smaller statistical errors, we test Casimir scal- 
ing through these. The most straightforward test is to 

note that (i^) must be independent of D if Casimir 
scaling holds. In Figure [7] we show that deviations from 
Casimir scaling are visible only very close to Tc. This 
has implications for weak coupling expansions. Beyond 
two-loop order in a perturbative series of the Polyakov 
loop, Casimir scaling violations can appear 40J|. These 
must be strongly suppressed compared to contributions 
which scale with the quadratic Casimir operator. 
An equivalent test is to note that the quantities Z]j{g^) 
which are obtained using the assumption in eq. (|14p are 
equal within errors, as shown in Figure [H Note that this 
agreement is an outcome of the renormalization proce- 
dure, and not built into it. A finer test of Casimir 
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FIG. 9: Difference loops for the sextet, 5l(,, and adjoint. Sis 
Polyakov loops using Casimir scaling psp . 
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scaling is obtained using the difference loops 

6Ld =L3- {Lof^^ , (15) 

The results for the renormalized difference loops for 
D = Q and 8 are shown in Figure [9] They are consistent 
with zero at all temperatures. For the higher representa- 
tions a statistically finer test is obtained with bare loops, 
since the errors on the renormalized loops are large. The 
results are shown in Figure 1101 Even here, Casimir scal- 
ing is a good approximation, which gets better the higher 
the temperature. 

A. Two-flavor QCD 




FIG. 10: Difference loops for all representation, &Id, using 
Casimir scaling (|15p . 



FIG. 11; Testing Casimir scaling for the bare Polyakov loop 
in 2-flavor QCD {Ld) from a 16^ x 4 lattice for all D = 
3,6,8,10,15. 

It is interesting to check whether Casimir scaling also 
holds for QCD with quarks. Since dynamical quarks 
break the center symmetry explicitly, and Nf — 2 
QCD has a finite temperature cross over rather than a 
true phase transition, the thermodynamic limit of the 
Polyakov loop below Tc is non-vanishing. 

In Figure [m we show {L^of'^" for D = 3,6,8 at all 
temperatures, and D = 10, 15 above Tc- The latter rep- 
resentations are too noisy below Tc to add any informa- 
tion. These scaled quantities are almost independent of 
D down to ^ l.bTc- Below this temperature we observe 
deviations to smaller values for the fundamental repre- 
sentation, whereas the values for higher D still coincide 
within errors. Therefore we see a violation of Casimir 
scaling, between the fundamental and other representa- 
tions, when entering the transition region which contin- 
ues to the smallest temperatures analyzed. These viola- 
tions are relatively mild. Differences between (L^)^/''° 
for D — 6 and D = 8 remain statistically insignificant 
even at the smallest temperatures, as shown in the inset 
of Figure [TT] 
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THE LARGE-iV LIMIT 
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each spatial site on the lattice, x, contains a matrix val- 
ued "spin", l{x), 



Z = J Y\_dl{x)exp[-SMM], 



(19) 



and the integration measure is the Haar measure. For 
SU{N) gauge groups the matrix takes values in SU(N). 
The action for these models can be written in the form 



Smm 



DD' 



0D,D'St{DxD'),O 

X Re^D (£')/£)/ {x + n) 

xh 

+lD6t{D),0^^elD{3 



(20) 



FIG. 



12: Differences loops for the sextet, 5/6, and adjoint. 



The relation between Polyakov loops in different rep- 
resentations becomes rather simple in the limit of a large 
numbers of colors, N [8]. In this limit the expectation 
value of a Polyakov loop in representation D can be ex- 
pressed in powers of the fundamental (Ln) and anti- 
fundamental, {Lj^ ~ (Ln)*) loop. 



tP+ tP- 



o 



(16) 



where the integers p+ and p_ are determined from the 
Young tableaux of the representation D. We investigate 
this large-iV factorization using our data obtained for 
TV = 3. 

Following we analyze difference loops 



(17) 
(18) 



Naively, the correction terms are expected to be of the or- 
der of (^3)^/3, i.e., about 33%. Results, shown in Figure 
[13 are clearly non-zero, except at around 2.5 where all 
loops are one. Our results are comparable in magnitude 
to those in but show rather different temperature de- 
pendence. Below 2Tc the corrections are relatively large, 
and the usefulness of the large-A^ approximation seems 
doubtful. However, above this temperature, the differ- 
ence loops are of order 5-10% of the loop itself, and 
therefore significantly smaller than the naive expecta- 
tions. The large-iV approximation seems to fare better 
than expected. This is similar to the conclusion reached 
for the equation of state in Q. 



VI. MATRIX MODELS 

One could seek effective field theories for Polyakov 
loops in the form of matrix models, i.e., models in which 



where x runs over every site in the lattice, n over the 
2d nearest neighbors. Id is the Polyakov loop in the ir- 
reducible representation D, i.e. the trace of the matrix, 
and t{D) is the triality of the irreducible representation 
D. In this section we use the notation of Appendix [A} 
i.e. traces are normalized to the dimension of the cor- 
responding representation, D. We also use the notation 
— Hd)- The constraint of vanishing triality arises 
from the Zn center-invariance of the pure gauge theory. 

For SU(3), the effective action with only the lead- 
ing term /Sg 3 has been investigated extensively over the 
years. However, when adding all irreducible representa- 
tions upto a certain £), as D varies, one needs the cou- 
plings 
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03,3 
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03,6, 0()fi 
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08,S, 78 


10 


08,10, /?10,10: 710 


15 


03,15' 0e.i5 


15' 


03,15''' 0S,15', 015 J5i' 



(21) 



and so on. 

A matrix model would be used to obtain the loop expec- 
tation values 

aB({/3})> = | Jl[dl{x)lDeM-SMM{l,{0})], (22) 



where {/?} denotes the whole set of couplings in the ac- 
tion. Equating these expressions to a sufficient number 
of observations on ^^(r), one would obtain the temper- 
ature dependence of the couplings. Other predictions of 
matrix models, which we do not explore here, are expec- 
tation values of moments (for example, the Polyakov loop 
susceptibilities) and correlation functions of loops. 
Note an intrinsic complication in the matching proce- 
dure. Since ^^(T) is scheme dependent, the couplings 
that one extracts by any matching procedure must also 
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values can be traded for the couplings, and all other ex- 
pectation values can be written in terms of these. For 
example, for a one-coupling matrix model, Cjv = 1, one 
could write ^^(^3), for all D > 3. This relation is RG in- 
variant: in two different renormalization schemes, if the 
values of ^3 at two different temperatures are the same, 
then the values of will be equal, for each D. For the 
pure gauge theory the data in Figure [T3] shows that the 
SU(3) pure gauge theory requires a matrix model with 
only a single coupling. A single coupling matrix model 
means that ratios such as /33,6//38,8 a-re fixed, and only 
one coupling is dependent on the temperature. 
The temperature independent ratios of couplings define 
the shape of the universal curves, ^^(^3), and the sin- 
gle tunable coupling says how the curve is traversed, in 
a given renormalization scheme, as T changes. There- 
fore, one can solve the problem in two steps: first use 
the universal curves to fix the ratios of the couplings, 
and finally solve the easier problem of finding the single 
left over coupling. Since exact solutions for the loop ex- 
pectation values are not known for matrix models with 
Nc = 3, one has to either solve the problem through a 
Monte Carlo simulation or in mean field theory. Here we 
investigate the latter option. 

Taking into account the irreducible representations 3 and 
6 in the effective action, one has 



s 



-3E 



FIG. 13: Different renormalized Polyakov loops shown as a 
function of the fundamental loop. The fact that the data 
in several different renormalization schemes (^3(12Tc) — 0.5 
in red, 0.75 in green and 1 in blue) collapse on to a single 
universal curve in each case implies that there can be only a 
single coupling matrix model which describes this data. The 
line is the result of a fit using four terms in the action, as 
described in the text, tuned to bracket the observed curve 
4(4). 



be scheme dependent. Note also that, in order to make 
contact with a matrix model, one has to choose a renor- 
malization scheme in which Pjj{T) < D. To the best 
of our knowledge, these points have not been noted in 
the literature. Furthermore, as the number of irreducible 
representations increases, the number of couplings in the 
effective theory which need to be matched to data in- 
creases rapidly. A determination of the effective action 
involves extraction of the couplings through such match- 
ing at each temperature. This is an ill-conditioned prob- 
lem unless the series can be cut off, and the number of 
couplings required is less than the number of pieces of 
data. 

It is interesting to ask how one can bound the number of 
couplings needed in the matrix model. If there are Cjv 
couplings to be determined, then Cn oi the expectation 



+p3,eRel3{x)le{x + n) 



+pQ-^Rele{x)re{x + h) 



(23) 



where Iq = (^3)^ — l^. Using this SU(3) relation, and 
making a mean-field approximation, we find that 



S = -6dV 



P^^^esReh + ^{eeRek + hReHl - l^)} 



+l3,^e,Re{ll - i;) 



(24) 



It is also possible to extend such a mean field treatment 
to models which include the octet representations. Using 
the invariance of the Haar measure, we can diagonalize 
the matrix, so that ^3 = exp(i(/))-t-exp(i7/))-|-exp[i((/>-|-'0)], 
and perform the integration over the remaining variables 
to give 

dh - -i^{l-cos(0- V')}{l-cos(20-KV)} 

x{l-cos(0 + 2V')}d(/)d'0- (25) 
Putting all this together, we find 

J dhexp{-S/V) 
exp[-UF]. (26) 



^(/?3,3>/33,6,/36,6>^3,4) 



1 V 
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In the mean-field theory the expectation values are com- 
puted simply as 



(27) 



where the Id can be expressed in terms of the angles ip 
and (j> using the formulae in eqs. (|A12IIA18|) and the re- 
lation Re/3 = cos ip + cos -I- cos(V' — (j)). 
Some of the results are shown in Figure [131 The observed 
Casimir scaling of loops implies a power-law dependence 
of loops on each other. The matrix model which includes 
only the coupling (3^ 3 is in fair agreement with the uni- 
versal curve £^{£3) when £3 > 1/2. However, it disagrees 
with the curve for Isi^s) already when £3 — 0.9. By in- 
cluding terms in /Sa^g and /3g g (in a fixed T-independent 
ratio to 3) the curve for £q can be improved; but this 
leads to no perceptible change in the curve for Is- How- 
ever, by introducing the coupling /3g g, and tuning the T- 
independent ratio P^-g/ ^■^3, one can contrive to improve 
the description of the two universal curves. However the 
universal curves for £iq etc., need further tuning. 
The conclusion seems robust: the SU(3) pure gauge the- 
ory data can be described within a single coupling matrix 
model. However, it seems hard to construct a matrix 
model with a small number of terms which reproduces 
the power-law dependence of £d on £3. The second re- 
sult has been obtained within a mean field theory, and 
needs verification in a more complete approach, such as 
the full simulation of such matrix models. 



VII. ADJOINT SOURCES AND GLUELUMPS 

Polyakov loops in representations with non-zero trial- 
ity vanish in the confined phase of the pure gauge theory, 
since the Z(3) symmetry of the action is realized on the 
states with large contribution to the path integral. The 
behavior of loops with vanishing triality can be different, 
because they are blind to the Z(3) symmetry involved in 
the QCD phase transition. This study is confined to the 
octet loop since all other triality zero loops that we con- 
structed turned out to have very large errors below Tc- 

A dynamical picture has been advanced for the behav- 
ior of adjoint Polyakov loops below Tc- An adjoint source 
can couple to a gluon in the medium to form a colorless 
composite called a gluelump. Correlations of triality zero 
loops can clearly be mediated by gluon exchange at any 
temperature, leading to screening. Gluelumps provide 
a summary of the main features of such screening j45| 
through two parameters: the free energy of separated 
gluelumps determines the asymptotic value of an octet 
Polyakov loop, the string breaking distance is the dis- 
tance at which the long-distance screening behavior sets 
in. 

We found that the bare adjoint Polyakov loop L\ below 
Tc assumes its thermodynamic limit, i.e., becomes inde- 
pendent of the volume for x 4 lattices with > 24. 
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FIG. 14: The renormalized adjoint Polyakov loops below and 
above the critical temperature compared to the fundamental 
loops above Tc. The solid line below Tc indicates the vanishing 
fundamental loops in the confined phase. 



Lg could be renormalized using either the QQ method for 
the octet loop explained in Appendix|B] However, we cal- 
culated L| at more couplings than the adjoint correlator. 
Since we found that agrees with Z3 (see Appendix iBjl . 
we used the .^3(17^) given in Table HI] to obtain L^(T). 
Table HTll lists the values found for the renormalized ad- 
joint Polyakov loop Lg and Figure [HI compares these re- 
sults to those for the renormalized fundamental Polyakov 
loop, L3. Although Lg becomes rather small below Tc, it 
is clearly non-zero for all temperatures analyzed by us. 
We observe that rises from 0.0087(16) at T/Tc = 0.907 
to 0.0219(48) just below T^ at T/Tc = 0.995. Cross- 
ing the critical temperature into the deconfined phase, 
Lg jumps almost an order of magnitude to 0.154(37) at 
T/Tc = 1.005. It is a little surprising to find the octet 
loop, which is blind to the Z(3) symmetry, change discon- 
tinuously at the symmetry-breaking transition. However, 
other triality-zero operators also change discontinuously 
at the phase transition, most notably the energy density. 
We now address the issue of string breaking and deter- 
mine the binding energy of the gluelump. The specific 
situation at T/Tc = 0.959, shown in Figure [T51 serves as 
an example. We have shown F^q g, and, since it becomes 
too noisy at large distances, also the "color average" free 
energy Lqq g, which has the same value at long distances. 
These free energies clearly show that adjoint sources are 
screened at large distances, in clear contrast to the lin- 
early rising free energy, Tqq 3 of sources in the funda- 
mental representation. We calculate the free energy at 
infinite separation between the sources using the cluster 
property, 

i^8,oo(T) = -2TlnLg(T) = 2mgiu,„p(T). (28) 

At zero temperature the energy stored in the field suf- 
fices to put on-shell two gluons from the medium and 
form two disjoint gluelumps ^4S\. At finite temperature 
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FIG. 15: Heavy quark-antiquark free energies for adjoint 
sources in the color singlet and color averaged channel com- 
pared to the Casimir scaled color singlet free energy of fun- 
damental sources at a temperature of 0.959 Tc. The dashed 
line indicates the definition of the string-breaking radius as 
explained in the text and the horizontal solid lines show the 
asymptotic value for the adjoint free energies. The thick black 
line indicates the adjoint T — potential Vs{r) = dsV3{r). 

this free energy can be identified with twice the gluclump 
screening mass. 

Results for Fs^cx. are collected in IIIII Fg.oo, is shown in 
FigurelTBl (upper panel). It changes httle with T, starting 
from 2.331(88) GeV at T/Tc = 0.907 and subsequently 
faUing to 2.06(12) GeV just below T^. At the lowest 
temperatures discussed here, Foo indeed approaches the 
T = value of Vg.oo = 2.4 - 3.0 GeY, which is twice the 
mass of the gluelump obtained in [43| . 
We define the string breaking distance rgtring (T) by com- 
paring the Casimir scaled free energy with fundamental 
sources, ^sF^— ^, with the screened value of the free en- 
ergy with adjoint sources, 

daF^Q^r^tringiT)) = FsMT). (29) 

The results are collected in Table lllll In FigurefTBl (lower 
panel) we show the resulting values of rgtring as a function 
of the temperature. There is rather mild change in rgtring 
with T. It varies from 1.180(61) fm at T/Tc = 0.907 to 
1.053(81) fm just below Tc- At the smallest temperature 
''string almost coincidcs with the T = value of 1.2 fm 

VIII. CONCLUSIONS 

We examined the renormalized Polyakov loop in many 
different irreducible representations of the gauge group 
SU(3) in the thermodynamic limit of pure gauge QCD. 
It has been known for a long time that the ultraviolet 
divergences of the Polyakov loop can be absorbed into a 
multiplicative renormalization "constant" Z(g^), where 



FIG. 16: Asymptotic values of the adjoint heavy quark free 
energies (upper panel). Estimates of the string-breaking ra- 
dius using the two definitions explained in the text (lower 
panel) . 





Foo [GeV] 


rstring(V8) [fm] 


LI 


0.907 


2.331(88) 


1.180(61) 


0.0087(16) 


0.924 


2.310(84) 


1.170(58) 


0.0099(17) 


0.942 


2.274(76) 


1.153(53) 


0.0116(17) 


0.959 


2.204(70) 


1.121(48) 


0.0143(19) 


0.977 


2.161(73) 


1.101(50) 


0.0168(23) 


0.986 


2.09(14) 


1.069(97) 


0.0198(53) 


0.995 


2.06(12) 


1.053(81) 


0.0219(48) 


0.000 


2.4-3.0 


~ 1.2 





TABLE III: Temperature dependence of Foo, string breaking 
distance ratring for the adjoint singlet free energy with respect 
to Vs (see text) and the renormalized adjoint Polyakov loop 
I/g. The last line gives the values at T = for twice the mass 
of the gluelump [47| and for the string breaking distance [isl ]. 



is the bare coupling. Such a renormalization fac- 
tor does not depend on long distance physics, such as 
the temperature, T (see eq. [7]). We implemented such a 
renormalization procedure by explicitly constructing an 
iteration using only explicitly gauge invariant quantities 
starting from a reference temperature T^ef incorporating 
this idea (see Section UlIB I for details). This so-called di- 
rect renormalization procedure was then used to extract 
the renormalized Polyakov loops in representations upto 
the 27 of SU(3) for a wide range of temperatures (see 
Section ini)- 

The technical part of our paper also consists of extend- 
ing the QQ renormalization procedure of [l^j to Polyakov 
loops in arbitrary representations of the gauge group (see 
Appendix [B|) . This is done by matching (gauge variant) 
correlation functions of sources in arbitrary representa- 
tions to zero temperature values at the ultraviolet cut- 
off. Although one does not demand explicitly that the 
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renormalization constant depends only on the bare cou- 
pling, the matching to zero temperature in the ultraviolet 
makes sure that this occurs. We checked that both renor- 
malization procedures have the same, one real parameter, 
freedom of choice of scheme (see Appendix [C|) . Having 
two drastically different renormalization schemes which 
are functionally equivalent allows us not only to use the 
most convenient scheme in any situation, but also to cross 
check the results by using both schemes whenever possi- 
ble. This puts the results of the lattice measurements on 
very strong footing. Furthermore, the equivalence of the 
two procedures shows that the short distance as well as 
the large distance parts of the heavy quark free energies 
obtained in Coulomb gauge become gauge independent 
as proposed in [igj. 

An interesting simplification occurs when Polyakov 
loop expectation values satisfy a relation called Casimir 
scaling (eq. E]) . Then the renormalization factors in all 
the different representations essentially boil down to a 
single factor. Furthermore, large- iV factorization evolves 
from the large-iV limit of the quadratic Casimirs, and 
hence Casimir scaling could provide an alternative route 
to large- iV scaling. We have presented tests of Casimir 
scaling in Section IIVI and of direct large- factorization 
in Section |Vl Both turn out to be reasonably reliable 
away from Tc. However, Casimir scaling is significantly 
more reliable and may provide a good route to scale large- 
A^ predictions down to A^ = 3. 

One subject of abiding interest is whether the high 
temperature phase of QCD can be described by a matrix 
model. We test this question in Section |Vll Casimir 
scaling implies that there are universal (renormalization 
scheme independent) relations between the renormalizcd 
Polyakov loop expectation values such that all the loops 
we studied depend only on the value of the fundamental 
loop. This implies that a matrix model description could 
work well away from T^. A single parameter variation 
of all couplings in the model would then reproduce the 
data on Polyakov loops, the temperature dependence of 
the couplings being, of course, renormalization scheme 



dependent. However, it seems that a simple model with 
a small number of parameters is not able to reproduce 
the power-laws in the lattice data, at least within the 
mean-field analysis of the matrix model performed here. 

Due to the Z(3)-symmetry of the pure gauge theory, 
all Polyakov loops with non-zero triality vanish in the 
confined phase of the pure gauge theory. For the ad- 
joint representation we have observed small, but non- 
zero, values below Tc for the first time in the thermody- 
namic limit (see Section|Vll|. Since static adjoint sources 
can form bound states, called gluelumps, with dynami- 
cal gluons, correlations of adjoint loops show screening 
(string breaking) even in the confined phase. As a re- 
sult, heavy quark free energies have a finite asymptotic 
value while for zero-zero triality they rise linearly with 
distance. Some aspects of the free energy can be cap- 
tured into the phenomenology of gluelumps through a 
mass and radius parameter. We present results for these 
quantities. 

Our primary technical result is the systematic devel- 
opment of two parallel renormalization procedures for 
Polyakov loops in arbitrary representations of the gauge 
group. This allows us to check that Casimir scaling of 
the rcnormalized loops is satisfied to good accuracy away 
from Tc. This is our main physical result, since it leads on 
to the discussion of large- A^ factorization and the matrix 
model description of lattice data. 
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APPENDIX A: POLYAKOV LOOPS IN 
IRREDUCIBLE REPRESENTATIONS OF SU(3) 

In order to obtain Polyakov loops in higher irreducible 
representations of SU(3) than the fundamental, we may 
use the theorem that the character in a direct product 
representation is the product of the corresponding char- 
acters, XPxgig) = Xp{9)Xq{9)- Then the direct product 
can be reduced using the Clebsh-Gordan series to yield 
the Polyakov loop in various representations. 

We use the following identities: 



(Al) 
(A2) 
(A3) 
(A4) 

+ (1,0) (A5) 
(A6) 
(AT) 

+ (0,2) (A8) 

+ (0,0) (A9) 

where we have specified the irreducible representations 
both in terms of its dimension and through the canon- 
ical label {p, q) where p and q are integers. Recall that 
the maximum weight in irreducible representation {p, q) 
is m = {(j3+g)/2V3, {p—q)/6}, and the dimension of this 
irreducible representation is D = (p+l)(g+l)(p+g+2)/2. 
The notation 15' stands for the irreducible representa- 
tion (4,0), and 15 denotes the irreducible representa- 
tion (2,1). Note that interchanging p and q gives the 
complex conjugate irreducible representation. The trial- 
ity of an irreducible representation can be defined to be 
t = (p ^ q)\:i. In each expression above, the trialities of 
all the irreducible representations on the right must be 
equal, and must equal the sum of the trialities (mod 3) 
of the irreducible representations on the left. This can 
be used as a check. 

More concretely, take the product over links and the trace 
defining the Polyakov loop in the 3 of SU(3), 





3x3 


= 


6 + 3 


(1,0) 


X (1,0) 


= 


(2,0) + (0,1) 




3x3 




8 + 1 


(1,0) 


X (0,1) 


= 


(1,1) + (0,0) 




6x3 


— 


10 + 8 


(2,0) 


X (1,0) 




(3,0) + (1,1) 




6x3 




15 + 3 


(2,0) 


X (0,1) 




(2,1) + (1,0) 




8x3 




15 + 6 + 3 


(1,1) 


X (1,0) 




(2,1) + (0,2) 




10 X 3 




15' + 15 


(3,0) 


X (1,0) 




(4,0) + (2,1) 




10 X 3 




24 + 6 


(3,0) 


X (0,1) 




(3,1) + (2,0) 




6x6 




15' + 15 + 6 


(2,0) 


X (2,0) 




(4,0) + (2,1) 




6x6 




27 + 8 + 1 


(2,0) 


X (0,2) 




(2, 2) + (1,1) 



hix) 



Tr l[Utia 

n=l 



ni), 



(AlO) 



where ^3(2;) is a complex number. Here the trace is 
normalized such that the unit matrix traces to 3. The 
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D 


(p. a) 


t C2{D) 


do 




3 


(1,0) 


1 


4/3 


1 




3 


(0,1) 


2 


4/3 


1 




6 


(2,0) 


2 


10/3 


5/2 




8 


(1,1) 





3 


9/4 


Im(L8) = 


10 


(3,0) 





6 


9/2 




15 


(2,1) 


1 


16/3 


4 




15' 


(4,0) 


1 


28/3 


7 




24 


(3,1) 


2 


25/3 


25/4 




27 


(2,2) 





8 


6 


Im(L27) = 



TABLE IV: Quadratic Casimir C2{D) for the representation 
D of SU(3), t — p — q mod 3 is the triality. do is the ratio 
C2{D)/C2{3). 



Polyakov loop in an irreducible representation is the char- 
acter in that irreducible representation. Hence, given the 
loop in one irreducible representation, that in the com- 
plex conjugate irreducible representation is obtained by 
complex conjugation. One specific example is 

%(x) = lUx), (All) 

where I* is the complex conjugate of L The Polyakov loop 
in the trivial irreducible representation 1 is unity (which 
gives vanishing potential in this irreducible representa- 
tion). Next we construct the series of other Polyakov 
loops. 



leix) 




(A12) 


hix) 


= Mx)f-i, 


(A13) 


hoix) 


= l3{x)l6{x) - Isix), 


(A14) 


hsix) 


= ll{x)l6{x) - laix), 


(A15) 


ll5' [x) 


= l3ix)ho{x) - hsix), 


(A16) 


h^ix) 


= llix)ho{x) - leix), 


(A17) 


hrix) 


= \leix)\^ - Isix) - 1. 


(A18) 



Further irreducible representations can be obtained if 
needed. Two of the reductions for the direct product 
have not been used. The Polyakov loop values in nor- 
malization used elsewhere in this paper is obtained by 
writing Luix) = lDix)/D. 



APPENDIX B: QQ RENORMALIZATION 

The QQ renormalization procedure [23 | can be ex- 
tended to static sources in arbitrary representations of 
the color group. For simplicity we will only discuss ad- 
joint sources here in detail, but the generalization to 
other representations is straight forward. 
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FIG. 17: Comparison of the color singlet quark-antiquark 
free energies for fundamental and adjoint sources scaled by 
the corresponding Casimir factor. The solid line represents 
the zero temperature potential, Vsir) — V3(r)/C2(8). 
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FIG. 18: The renormalization constants, Zoig^), for funda- 
mental {D — 3) and adjoint {D = 8) sources plotted vs. the 
bare coupling — 2N/f3. The solid line is the same as in 

fig.H 

Given an SU(3) matrix in the fundamental representa- 
tion, , the corresponding adjoint matrix, is 

C/8 =iTr(A,;73A,C/-^^), i,j^l,...,8, (Bl) 

where are the Gell-Mann matrices. From the hermitic- 
ity of Xi and cyclicity of the trace all matrix elements Ufj 
are real. This formula can be used to convert all link el- 
ements from the fundamental to the adjoint. 

The thermal Wilson line in the adjoint representation 

is 

Psix) = Y[ U!ix,xi). (B2) 

X4 — Q 

Another way to define it is to take the fundamental Wil- 
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son line of eq. ((T|) and convert it to the adjoint using the 
prescription of eq. (|Bip . The adjoint Polyakov loop is 
the trace 



(B3) 



As before, we have normalized the trace such that the 
trace of the unit matrix is 1. 

Define the correlator of two adjoint thermal Wilson 
lines, 



(B4) 



and r ~ \xi — X2\- This correlator is clearly gauge de- 
pendent, and hence we define it through Coulomb gauge 
fixing (see [i^lHO] for more on this point). The free en- 
ergy with two static adjoint sources a distance r apart 
is 



(B5) 



Here C and F denote bare correlators and free energies; 
the same notation without a tilde will denote renormal- 
ized quantities. 

Since the Polyakov loop is renormalized multiplica- 
tively, the free energies are additively renormalizable. We 
match the finite temperature free energy to the zero tem- 
perature potential at the smallest attainable distance, a, 
on a lattice 



QQ,8 



a, T) + 2Tds In 



. 

where Vsir) is the zero temperature potential in the 
fundamental representation. We use the potential de- 
rived in [sij . In the matching procedure we have used 
Casimir scaling of the potential at short distances. This 
is seen in continuum [40| and lattice [4l| perturbation 
theory. In Figure [17] we show the renormalized quark- 
antiquark free energies for static sources in the adjoint 
{D — 8) and fundamental {D — 3) representations to- 
gether with the Casimir scaled zero temperature poten- 
tial Vs{r) = dsV3{r). The data clearly validates the as- 
sumption of short distance Casimir scaling on which the 
procedure rests. 

Once the free energies have been renormalized at small 
distances, their large distance behavior is fixed. The 
asymptotic value can be used to the define the renor- 
malized Polyakov loop through the cluster property 



lim JCl^^Jr,T) 



lim exp 



2T 



(B7) 



This completes the QQ renormalization procedure for the 
octet loop. The requirement that Zg, depend only on 
is not explicitly imposed in the QQ renormalization pro- 
cedure. However, the results, plotted in Figure [T51 show 



that this is obtained. The figure also shows that Casimir 
scaling of the renormalized Polyakov loop is obtained, 
since Zs and Z3 agree. 

Any 3x3 unitary matrix with unit determinant is 
uniquely specified by eight real numbers, which are coor- 
dinates in the abstract group SU(3). Given these coor- 
dinates, there are canonical techniques for building ma- 
trices in any representation D which generalize cq. (|B1|) . 
Hence, given the thermal Wilson line in the fundamental, 
one can construct the equivalent for arbitrary D. From 
that one can generalize every step of the procedure from 
eq. (jB4[) on for any D. 



APPENDIX C: RENORMALIZATION SCHEMES 
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FIG. 19: Renormalization constants for fundamental loops 
at different scales C = 0.0, —^/a (upper data points) and 
^/a (lower data points) in (IC1|I from the direct and QQ- 
renormalization method. The solid lines are the properly 
scaled fit from fig. [l] 



In the direct renormalization procedure the free- 
dom of scheme choice is the multiplicative ambiguity 
L^nC^ref) — > KjjL'^jjiTref), for some constant Ku. This 
implies that at another temperature the renormalized 
Polyakov loop is scaled by the factor K^"^^'^ . In the 
QQ renormalization procedure it is the freedom of defin- 
ing the zero of the T = potential 



Vz{r) ^Vz{r) + C. 



(CI) 



Using Casimir scaling for the short distance potential, 
this clearly leads to the scaling 



Ld{T) 



(C2) 



which incorporates Casimir scaling for the Polyakov loop. 
Our standard scheme choice corresponds to choosing C 
such that the triplet potential at T = is given by the 
results of I5ll. 
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In Figure [19] we show examples of the change in renor- 
mahzation scheme using the two procedures. Despite the 
scaUng freedom, the dependence of Z{g'^) on the bare 
couphng is independent of this scale and the tempera- 
ture dependence of the renormalized Polyakov loops fol- 
low from El. 



APPENDIX D: COLOR AVERAGE FREE 
ENERGIES 
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FIG. 20: Average free energies for D — 6 (top) and D — 8 
(bottom) compared to the fundamental free energy below Tc. 

We now turn to color averaged QQ-free energies. Since 
the corresponding correlators can be obtained without a 
costly gauge fixing, we were able to calculate FQQ jj{r, T) 
for representations D = 3,6,8 in the temperature range 
0.9 — 3 Tc on 32^ x 4 lattices. The color average correlator 
for at temperature T in representation D is defined by 



CQQ,D{r,T) = {L'-j,{x,W^{x^))), 



(Dl) 



where the star denotes complex conjugation. The color 
average free energy FQQ jj{r,T) = -T log Cqq_^ (r, T). 
If Casimir scaling holds, then 

FQQ,D{r,T)IC2{D) = FQQ,j,,{r,T)IC2{D'). (D2) 
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FIG. 21: Average free energies for D = 6 (top) and D = 8 
(bottom) compared to the fundamental free energy above Tc. 



We test this relation below and above Tc- 

Below T(. we employ the renormalized average free en- 
ergies 

FQQ^D{r,T)^FQQ^u{r,T) + 2TdD\nZD{g^) (D3) 

for the three lowest temperatures divided by their 
Casimir in fig. 1201 The renormalization constants used 
are those found by the renormalization procedure de- 
scribed in sec. IIII CI Thus for the smallest distances all 
curves coincide as a consequence of the renormalization 
procedure. However, for all T < Tc and representations 
D = 6,8 deviations to smaller values start to show up 
quite early, i. e. for separations ry/a>l for D — 6 and 
rY^>0.8 for D — 8, respectively. This effect is more pro- 
nounced for the adjoint average free energy than for the 
sextet. The effect of string breaking sets in at larger dis- 
tances than shown here and will be discussed in sec. IVIII 
Above Tc we compare 



AF 



(r, T) = FoQ^jjir, T) - FooAr -> ^, T) (D4) 



divided by their Casimir for the same representations 
D = 3, 6, 8 in fig. [51] We observe screening to take 
place in both higher representations. The curves for both 
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D = 6 and D = 8 deviate to smaller values compared to 
the fundamental case. We find that the ordering 

^FqqA^^T) AFQQ^,{r,T) AFQQ^,{r,T) 
C2(6) C2(8) ^ 6-2(3) 

(D5) 



holds throughout the entire distance interval above Tc- 

Thus, we conclude, that Casimir scaling (jD2p is clearly 
violated for the average QQ free energies in the temper- 
ature range 0.9 — 3Tc for the fundamental, sextet and 
adjoint representations. 



